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Using the QGP motivated threshold model, where all the J/tp's are suppressed above a threshold 
density, we have analyzed the preliminary PHENIX data on the centrality dependence of nuclear 
modification factor for J/tp's in Cu+Cu and in Au+Au collisions, at RHIC energy, ^^^=200 GeV. 
Centrality dependence of J/tp suppression in Au+Au collisions are well explained in the model for 
threshold densities in ranges of 3.6-3.7 fm~ 2 . J/tp suppression in Cu+Cu collisions on the other 
hand are not explained in the model. 
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PACS numbers: PACS numbers: 25.75.-q, 25.75.Dw 



I. INTRODUCTION 

In relativistic heavy ion collisions J / tp suppression has 
been recognized as an important tool to identify the pos- 
sible phase transition to quark-gluon plasma. Because 
of the large mass of the charm quarks, cc pairs are pro- 
duced on a short time scale. Their tight binding also 
makes them immune to final state interactions. Their 
evolution probes the state of matter in the early stage 
of the collisions. Matsui and Satz p|, predicted that in 
presence of quark-gluon plasma (QGP), binding of a cc 
pair into a J/tp meson will be hindered, leadingto the so 
called J /if) suppression in heavy ion collisions [lj . Over 
the years, several groups have measured the J/tp yield 
in heavy ion collisions (for a review of the data prior to 
RHIC energy collisions, and the interpretations see Refs. 
0,0). In brief, experimental data do show suppression. 
However, this could be attributed to the conventional 
nuclear absorption, also present in pA collisions. 

In recent Au+Au collisions at RHIC, one observe 
a dramatic suppression of hadrons with high momen- 
tum, transverse to beam direction (high pt suppression) 
Q, y> H> Q- This has been interpreted as evidence for 
the creation of high density, color opaque medium of de- 
confined quarks and gluons @. It is expected that high 
density, color opaque medium will leave its imprint on 
J/tp production. At RHIC energy, it has been argued 
that rather than suppression, charmoniums will be en- 
hanced 0, IToj . Due to large initial energy, large num- 
ber of cc pairs will be produced in initial hard scatter- 
ings. Recombination of cc can occur enhancing the char- 
monium production. PHENIX collaboration have mea- 
sured the centrality dependence of J/ip invariant yield 
in Au+Au collisions at RHIC energy, y/s NN —200 GeV 
[Til IT^ . More recently they have improved upon the 
statistics and preliminary results for the centrality de- 
pendence of nuclear modification factor (Raa) f° r J/ip 
in Au+Au collisions and in Cu+Cu collisions are avail- 
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able [T3. IT^j. In most central Au+Au collisions, J ftp's 
are suppressed by a factor of ~ 3. PHENIX data on 
J/ip production in Au+Au/Cu+Cu collisions, are not 
consistent with models which predict J/tp enhancement 
0, 0. It was also seen that various models, e.g. co- 
mover model 0| , statistical coalescence model or 
the kinetic model |17| . fail to explain the (preliminary) 
PHENIX data on the nuclear modification factor for J/tp 
in Cu+Cu and in Au+Au collisions. The data are also 
not explained in the Glauber model of normal nuclear ab- 
sorption [l9(. Recently, in a QCD based nuclear absorp- 
tion model , we have analyzed the preliminary PHENIX 
data on J/tp suppression in Cu+Cu and in Au+Au col- 
lisions!^. In the QCD based nuclear absorption model 
|2ll |22|. cc pair, during its passage through a nuclear 
medium, gain relative 4-square momentum. Some of the 
pairs gain enough to cross the open charm threshold and 
are lost. The model explained the PHENIX data on the 
centrality dependence of J/tp suppression in Cu+Cu col- 
lisions at RHIC but failed for Au+Au collisions. It was 
concluded that in Au+Au collisions, J/tp are suppressed 
in a medium, unlike that produced in SPS energy nu- 
clear collisions or at RHIC energy Cu+Cu collisions. If 
in Au+Au collisions, J/tp's are suppressed in a decon- 
fined matter, PHENIX data should be explained in a 
QGP motivated model, like the threshold model p3ll2^|. 
Blaizot et al [2j|, |24| , proposed the threshold model to 
explain the NA50 data on anomalous J/tp su ppre ssion 
in 158 AGeV Pb+Pb collisions at SPS energy [25| . To 
mimic the onset of deconfining phase transition above a 
critical energy density and subsequent melting of J/tp's, 
J / tp suppression was linked with the local energy density. 
If the energy density at the point where J/tp is formed, 
exceeds a critical value (e c ), J/tp's disappear. 

In the present paper, in the threshold model, we have 
analyzed the preliminary PHENIX data on the centrality 
dependence of J/tp suppression in Cu+Cu and in Au+Au 
collisions. As it will be shown below, while the centrality 
dependence of J/tp suppression in Au+Au collisions are 
well explained in the model, that in Cu+Cu collisions 
are not. We have also analyzed the PHENIX data on the 
centrality dependence of pr broadening of J / tp 0, ^| . 
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No definitive conclusions can be obtained from the pr 
broadening data . 

The plan of the paper is as follows: in section II, 
we briefly describe the threshold model. PHENIX data 
on the centrality dependence of J/tp suppression are 
analyzed in section III. In section IV, we analyze the 
PHENIX data on px broadening of J ftp. Summary and 
conclusions are drawn in section V. 



II. THRESHOLD MODEL 

The details of the threshold model could be found in 
[23L l24j . It is assumed that fate of a J/tp depend on 
the local energy density, which is proportional to partic- 
ipant density. If the energy density or equivalently, the 
participant density, exceeds a critical or threshold value, 
deconfmed matter is formed and all the J ftp's are com- 
pletely destroyed (anomalous suppression). This anoma- 
lous suppression is in addition to the "conventional nu- 
clear absorption" . Transverse expansion of the system is 
neglected. It is implicitly assumed that J ftp's are absor- 
ped before the transverse expansion sets in. 

In the threshold model, number of J/tp mesons, pro- 
duced in a AA collision, at impact parameter b can be 
written as, 



AA 



'(b) = a«* J d 2 sT^(s)T^(b-s) 

^ Sanom (b, s) , (1) 

where T e f * (b) is the effective nuclear thickness, 

/CO p DO 

dzp(b, z)exp(-a abs I dz/p(b,z/)), 
-OO J Z 

^abs being the J/V'-Nucleon absorption cross-section. 
£anorn(b,s) in EqQ] is the anomalous suppression factor 
introduced by Blaizot et al. [23, [24| ■ Assuming that all 
the J ftp's get suppressed above a threshold density (n c ), 
the anomalous suppression can be written as, 



»(b,s) = 8(n(b,s) - n c ) 



(3) 



where n c is the critical or the threshold density. n(b, s) 
is the local transverse density. At impact parameter b 
and at the transverse position s, local transverse density 
it can be obtained as, 



density n c . With J/?/>-nucleon absorption cross-section 
<jj/^n=6A mb, NA50 data are well explained in the 
model with n c =3.7 fm~ 2 . Better fit to the data is ob- 
tained if the theta function (Eq[3J) is smeared, at the 
expense of an additional parameter. Later experiments 
l 26] indicate that J/^-nucleon absorption cross-section 
is ~ 4 mb, rather than 6.4 mb. NA50 collaboration also 
revised their data [27| • The revised NA50 data were also 
analyzed in the threshold model (28|. With a a b s ~ 4 
mb, large smearing of the threshold density is required. 
Large smearing of threshold density, effectively excludes 
formation of deconfmed matter at SPS energy. 




FIG. 1: Transverse density in Cu+Cu (left panel) and in 
Au+Au (right panel) collisions, for various values of the im- 
pact parameter, b=0,2,4.... (from top to bottom). The origin 
is at a distance, d — 6/(1 + Ra/Rb) from the center of the 
nucleus A. 



III. J/tp SUPPRESSION IN CU+CU/AU+AU 
COLLISIONS 

In the threshold model, fate of a J/tp is determined on 
the local (transverse) density. If the local (transverse) 
density exceeds the threshold density, J/ip's are com- 
pletely destroyed. In Fig^ for a number of impact pa- 
rameters, the transverse density, n(b, s) in Cu+Cu and in 
Au+Au collisions are shown. We have used the Woods- 
Saxon form for the density, 



n(b,s) = T A {s)[l-exp{-a NN T B (s-b))} 

+T B (b-s)[l-exp{-cT NN T A (s))] (4) 

Blaizot et a/ P, H fitted the NA50 data [H on 
the transverse energy dependence of J/tp suppression in 
158 AGeV Pb+Pb collisions and obtain the threshold 



P( r ) = 7— fr m y [d 3 rp(r)=A (5) 

1 + exp((r - R)/a) J 

with i?=4.456 (5.415) fm and a=0.54 (0.535) fm, for Cu 
(Au) nuclei. 

For central collisions, maximum transverse density in 
Cu+Cu collisions is ~2.63 fm~ 2 , while that for Au+Au 
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collisions is ~ 4.32 /to -2 . Then if J ftp's are anomalously 
suppressed, say, above a threshold density, n c —3.7 fm~ 2 , 
J ftp suppression in Cu+Cu collisions will not be affected 
as the transverse density never exceeds the threshold den- 
sity. In Au+Au collisions, on the other hand, J ftp's will 
be anomalously suppressed. In Au+Au collisions also, 
only in collisions where local density n(b, s) exceeds the 
threshold density, J ftp's will be anomalously suppressed. 
In all other collisions, J ftp's will be absorped only due 
to J/^-nucleon interaction. Then if J ftp suppression is 
measured as a function of impact parameter or equiv- 
alently, as a function of centrality of collisions, sudden 
change of slope will be observed. 
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FIG. 2: Preliminary PHENIX data on the participant num- 
ber dependence of nuclear modification factor (Raa) for J ftp, 
in Au+Au collisions. Raa for J ftp in d+Au collisions are 
also shown. The dashed line is the Glauber model predic- 
tion for Raa with a a bs=3 mb. The solid lines (bottom to 
top) are the threshold model predictions for threshold den- 
sity, n c =3.6,3.7,3.8,3.9 and 4 /m -2 . 

Such a sudden change of slope in J ftp suppression is 
observed by the PHENIX collaboration in Au+Au colli- 
sions. In Fig[3 preliminary PHENIX data on the central- 
ity dependence of nuclear modification factor (Raa) for 
J ftp, in Au+Au collisions are shown [HEil- PHENIX 
collaboration has taken data in two ranges of rapidity in- 
tervals, (i) -0.35 < y < 0.35 and (ii)1.2 < y < 2.2. Both 
the data are shown. Even though data points are few, 
sudden change of slope of Raa around N part ~ 150 is 
evident. PHENIX data for Raa in d+Au collisions are 
also shown in the Figure. J ftp's are suppressed in d+Au 
collisions also. Glauber model analysis indicate that in 
d+Au collisions, J/i/'-micleon absorption cross-section is 
small, <j a bs ~ 1 — 3 mb 0. In Fig|2 prediction for Raa, 
in the Glauber model of normal nuclear absorption is 
shown as the dash-dot-dashed line. The prediction is ob- 
tained with a a bs =3 mb. It is interesting to note that the 



Glauber model of nuclear absorption, with er a {, s =3 mb, 
explains the suppression in peripheral and mid-central 
collisions. Only in very central collisions, the Glauber 
model of nuclear absorption model predict much less sup- 
pression than observed by the PHENIX collaboration. 
Anomalous suppression in Blaizot's threshold model can 
provide the additional suppression required in very cen- 
tral collisions. Due to paucity of data, we do not at- 
tempt to fit the PHENIX data and extract the threshold 
density, n c . Rather, for a number of threshold density, 
?i c =3.6,3.7,3.8,3.9 and 4.0 /to -2 , we obtain predictions 
for Raa in the threshold model. In Fig|3 solid lines are 
the predictions in the threshold model. With anomalous 
suppression, J ftp's are strongly suppressed in central col- 
lisions and it is evident that for threshold density in the 
ranges of 3.6-3.7 /to -2 , the threshold model describe the 
data adequately well. 
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FIG. 3: Preliminary PHENIX data on the centrality depen- 
dence of nuclear modification factor (Raa) for J ftp in Cu+Cu 
collisions. Raa for J ftp in d+Au collisions are also shown. 
The solid lines (top to bottom) are Glauber model predic- 
tions with a a bs= 1,3,5 and 7 mb. 

We now analyze the PHENIX data on the centrality 
dependence of nuclear modification factor, in Cu+Cu col- 
lisions. The data are shown in FigOU We have not distin- 
guished between the mid-rapidity and forward rapidity 
data. In contrast to Au+Au data, Raa in Cu+Cu colli- 
sions do not show any sudden change in slope. Threshold 
model, which induces sudden change in slope, is not war- 
ranted by the data. Indeed, if the threshold density is in 
the ranges of 3.6-3.7 /m -2 , anomalous suppression will 
be ineffective in Cu+Cu collisions. As shown in Fig^ 
maximum transverse density reached in central Cu+Cu 
collisions is ~ 2.63 /to -2 , much less than the threshold 
density 3.6-3.7 /to -2 . 

As mentioned earlier, J ftp production in d+Au col- 
lisions at RHIC, require a a bs ~ 1-3 mb. Glauber 
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model predictions with a a b s =3 mb, also explains the 
PHENIX data on the centrality dependence of J/ip sup- 
pression in peripheral and mid-central Au+Au collisions. 
Only in very central Au+Au collisions, data demand 
anomalous suppression. However, in Cu+Cu collisions, 
Glauber model predictions with crabs =1-3 mb fails to 
explain the centrality dependence of J /ip suppression. 
In Fig|3 Glauber model predictions for J /ip suppres- 
sion in Cu+Cu collisions, for various values of J/ip- 
nucleon absorption cross-sections, <7 {, a =l,3,5, and 7 mb 
are shown. For cr aos =l-3, the model could explain sup- 
pression only in very peripheral collisions (N part =10-20). 
In mid-central and very central collisions, the model with 
Cafes =1-3 mb, underpredict the suppression. If the J/ip- 
nucleon cross section is increased, while the model ex- 
plains the mid-central and very central collisions, it fails 
to explain peripheral collisions. The analysis indicate 
that, with a single value J/^-nucleon absorption cross- 
section u a bs, the Glauber model of nuclear absorption, 
can not explain the centrality dependence of J /ip sup- 
pression in Cu+Cu collisions. J/ip suppression in Cu+Cu 
collisions at RHIC energy is more complex than envis- 
aged in the Glauber model of nuclear absorption. Indeed, 
a more complex, QCD based nuclear absorption model, 
does explain the centrality dependence of J/ip suppres- 
sion in Cu+Cu collisions adequately well |20). 



pression. Successful description of PHENIX data on J/tp 
suppression in Au+Au collisions then strongly support 
formation of deconfined matter in central Au+Au colli- 
sions. However, we note that the model neglects some 
very important effects, e.g. (i) feed back from ip 1 and x 
states and (ii) transverse expansion. A considerable frac- 
tions of J/ip's are from decay of ip' and x states. That 
part is completely neglected here. Threshold density for 
anomalous suppression of higher states, ip' and x should 
be less than that for a J/ip. Then presently estimated 
threshold density n c =3.6-3.7 fm~ 2 represent an upper 
limit of the threshold density. At RHIC, model studies 
indicate that in the deconfined phase, the system un- 
dergoes significant transverse expansion [2j|. The local 
transverse density is a key ingredient to the Threshold 
model. In an expanding system, local transverse density 
will be diluted. J/tp's which are anomalously suppressed 
in a static system, may survive in an expanding system 
due to dilution. Then, the presently estimated critical 
density 3.6-3.7 fm~ 2 will again represent an upper limit 
of the threshold density. 

IV. CENTRALITY DEPENDENCE OF p T 
BROADENING IN CU+CU/ AU+AU 
COLLISIONS 
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FIG. 4: Collision number dependence of the ratio Nab / o" 9 at in 
Au+Au collisions. The dashed line is the ratio in the normal 
nuclear absorption model, with a a i, s =3 mb. The solid lines 
are Nab / &gN in the threshold model, with threshold den- 
sity n c =3.6,3.7,3.8,3.9 and 4.0 fm~ 2 (from bottom to top) 
respectively. 



It is well known that in pA and AA collisions, the sec- 
ondary hadrons generally show px broadening [30L l3l| . 
Pt broadening of J / ip in Cu+Cu and in Au+Au collisions 
at RHIC energy v^jvw - ^00 Ge V, h as been measured by 
the PHENIX collaboration They measured the 

collision number dependence of square of transverse mo- 
mentum for J /ip. It is interesting to compare the thresh- 
old model predictions with the PHENIX data. 

The natural basis for the px broadening is the initial 
state parton scatterings. For J/ip's, gluon fusion being 
the dominant mechanism for cc production, initial state 
scattering of the projectile/target gluons with the tar- 
get/projectile nucleons causes the intrinsic momentum 
broadening of the gluons, which is reflected in the px 
distribution of the resulting J/ip's. Parameterising the 
intrinsic transverse momentum of a gluon, inside a nu- 
cleon as, 



f(q T ) ~ exp(-q T / < q\ >) 



(6) 



momentum distribution of the resulting J /ip in NN col- 
lision is obtained by convoluting two such distributions, 



exp(-p T / < p T > J Jx) 



(7) 



In the threshold model, it is implicitly assumed that 
J/ip's are absorped in a deconfined matter. The critical 
energy density for deconfined matter formation is pro- 
portional to the threshold density. Melting of J/ ip due to 
color screening is mimicked by the sudden on set of sup- 



where < p T >jvjv =< 9t > + < It > - I n nucleus- 
nucleus collisions at impact parameter b, if before fu- 
sion, a gluon undergo random walk and suffer N num- 
ber of subcollisions, its square momentum will increase 
to qi, —y q T + NSq, 5q being the average broadening in 
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each subcollisions. Square momentum of J/tp then easily 
obtained as, 



< Pt >ab (b) =< Pt > J Jt +5oN AB (b) (8) 

where Nab (b) is the number of subcollisions suffered by 
the projectile and target gluons with the target and pro- 
jectile nucleons respectively. 

Average number of collisions Nab (b) can be obtained 
in a Glauber model |3lj . At impact parameter b, the po- 
sitions (s, z) and (b — s, z') specifies the formation point 
of cc in the two nuclei, with s in the transverse plane and 
z, z' along the beam axis. The number of collisions, prior 
to cc pair formation, can be written as, 



N(b,s,z,z')= <jgN I (1zapa(s,za) (9) 

J — oo 

+<r gN / dz B pB(b - s,z') 



where a g N is the gluon-nucleon cross-section. Above ex- 
pression should be averaged over all positions of cc for- 
mation with a weight given by the product of nuclear 
densities and survival probabilities S, 
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FIG. 5: J/ip mean square transverse momentum as a function 
of collision number, in mid- rapidity p+p, d+Au, Cu+Cu and 
Au+Au collisions are shown. The solid line is the best fit, in 
the threshold model, to the combined Cu+Cu and Au+Au 
data. 



NAB(b) — / d 2 s / dzpA{s,z) / dz'psib- s,z') x 

J J — oo J —oo 

S(b,s,z,z')N(b,s,z,z')/ / d 2 s / dzpA(s,z)x 

J J — oo 

dz pB(b — s, z)S(b, s, z, z) (10) 



Centrality dependence of the ratio Nab/^qn, 
in Au+Au collisions, for the threshold densities, 
ri c =3.6,3.7,3.8,3.9 and 4.0 fm~ 2 , are shown in Fig0] (the 
solid lines from bottom to top) . We also show the ratio in 
the normal nuclear absorption model (the dashed line). 
Nab I VgN increases with centrality, more central the col- 
lisions, the gluons suffer more number of collisions. In 
normal nuclear absorption model, Nab / c gn continues to 
increase with centrality (or collision number). However, 
rate of increase slows down at more central collisions. 
A different behavior is obtained in the threshold model. 
For a fixed threshold density n c , N abI^qu exactly corre- 
sponds to normal nuclear absorption model, till a colli- 
sion number N c . Beyond N c , Nab/vqN hardly changes. 
It is understood. Beyond a N c , transverse density ex- 
ceeds the threshold density and J ftp's are completely de- 
stroyed. As Nab is weighted by the anomalous suppres- 
sion, it hardly changes beyond that collision number. 

Pt broadening of J/ip's in AA collisions depends on 
two parameters, (i) < p T >^/^, the mean squared trans- 
verse momentum in NN collisions and (ii) the product of 
the gluon-nucleon cross-section and the average parton 
momentum broadening per collision, ct^jv^o- < Pt >nn 



is measured in RHIC energy p+p collisions, < p T >nn— 
4.2 ± 0.7 GeV 2 . As gluons are not free, the other param- 
eter, <j g N$o is essentially non-measurable. Its value can 
be obtained from experimental data on pt broadening 
of J/ip. In SPS energy S+U/Pb+Pb collisions <7 9 n5q is 
estimated as 0.442 ± 0.056 GeV 2 HI- a gN S at RHIC 
energy is of interest. 

PHENIX data on the centrality dependence of mean 
square transverse momentum < p T >, in Cu+Cu and in 
Au+Au collisions are shown Fig[S| < p\ > in in p+p 
and in d+Au collisions are also shown. Quality of data 
is poor, few data points with large error bars. Evidently 
data do not show any evidence of pt broadening. Within 
the errors, < p T > in Cu+Cu and in Au+Au collisions 
agree with that in NN collisions, pt broadening of J/tp's 
is minimum at RHIC. 

To find <J g N$o at RHIC energy, we fit the combined 
Cu+Cu and Au+Au data set (individual Cu+Cu or 
Au+Au data points are few). We fix < p T >mn at the 
measured central value, < p\ >nn — 4.2 GeV 2 , and 
vary a g isr5o. < p T > in Au+Au or in Cu+Cu show very 
little dependence on the threshold density. In FigO best 
fit is obtained with threshold density n c =3.6,3.7,3.8,3.9 
and 4.0 fm~ 2 are shown. They can not be distinguished. 
Best fit is obtained with a gN S = 0.31 + 0.48 GeV 2 . Due 
to poor quality of the data, the a g N5o is ill determined. 
Estimated error is larger than the central value. We con- 
clude that PHENIX data can not determine the a g Nfio 
at RHIC energy. 
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V. SUMMARY AND CONCLUSIONS 

To conclude, in the QGP motivated threshold model, 
we have analyzed the preliminary PHENIX data on the 
centrality dependence of J/ip suppression in Cu+Cu and 
in Au+Au collisions. In the threshold model, in addi- 
tion to the normal nuclear absorption, J ftp's are anoma- 
lously suppressed, such that, if the local transverse den- 
sity exceeds a threshold density n c , all the J ftp's are ab- 
sorped. The model predicts a sudden change of slope 
in J /ip suppression as a function of centrality. Prelimi- 
nary PHENIX data on the centrality (participant num- 
ber) dependence of J ftp suppression in Au+Au collisions 
at RHIC are well explained in the model. It also re- 
produces the sudden change of slope around participant 
number, N part = 150. Estimated threshold density ranges 
between 3.6-3.7 fm~ 2 . PHENIX preliminary data on the 
centrality dependence of J/ip suppression in Cu+Cu col- 



lisions do not show any sudden change of slope, charac- 
teristic of the threshold model. Local transverse density 
in most central Cu+Cu collision is only 2.63/m -2 , much 
less than the estimated threshold density in Au+Au col- 
lisions. Anomalous suppression is not effective in Cu+Cu 
collisions. However, J/tp production in Cu+Cu colli- 
sions is not explained the Glauber model of normal nu- 
clear absorption. While very peripheral collisions require 
Cabs ~T-3 mb, more central collisions require a a bs ~6-7 
mb. It appears that J/ip production in Cu+Cu collisions 
is more complex than envisaged in the simple Glauber 
model of nuclear absorption. We have also analyzed the 
PHENIX data on pr broadening. The quality of data is 
not good enough for a definitive conclusion. Apparently 
at RHIC energy, J/ip's donot show any px broadening. 

In conclusion, present analysis strongly support decon- 
fined matter formation in central Au+Au collisions at 
RHIC energy ^^=200 GeV. 
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